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Abstract 

In this paper we consider the problem of transmitting a continuous alphabet discrete-time source over an AWGN 
channel. We propose a constructive scheme based on a set of curves on the surface of a 2A'^-dimensional sphere. Our 
approach shows that the design of good codes for this communication problem relies on geometrical properties of 
spherical codes and projections of A'^-dimensional rectangular lattices. Theoretical comparisons with some previous 
works in terms of the mean square error as a function of the channel SNR as well as simulations are provided. 

I. Introduction 

The problem of designing good codes for continuous alphabet sources to be transmitted over a channel with 
power constraint can be viewed as the one of constructing curves in the Euclidean space of maximal length and 
such that its folds (or laps) are a good distance apart. When the channel noise is bellow a certain threshold, a bigger 
length essentially means a higher resolution when estimating the sent value. On the other hand, if the folds of the 
curve come too close, this threshold will be small and the mean square error (mse) will be dominated by larger 
errors as consequence of decoding to the wrong fold. Explicit constructions of curves for analog source-channel 
coding were presented, for example, in |9| and lITOll . 

In this work, we will consider spherical curves in the 2iV-dimensional Euclidean space. We develop an extension 
of the construction presented in ||9] to a set of curves on layers of flat tori. Our approach explores the geometrical 
properties of spherical codes and projections of A^-dimensional lattices. In the scheme presented here, homogeneity 
and low decoding complexity properties were preserved whereas the total length can be meaningfully increased. 
Theoretical results as well as simulations show that this approach outperforms previous ones in terms of the tradeoff 
between SNR and mean square error. 

This paper is organized as follows: in Sections II and III, we introduce some mathematical/information-theoretical 
tools used in our approach. In Section IV we present a scheme to design piecewise homogeneous curves on the 
Euclidean sphere and describe the encoding/decoding process. In Section V we derive an extended version of the 
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27 Lifting Construction {T\ suitable to our problem and present some examples and length comparisons with some 

28 previous constructions. In the last sections we present theoretical and experimental analysis of the proposed scheme. 

29 II. Communication Framework 

30 The communication system we considered here is illustrated below. The objective is to transmit an input real 

31 value X over a channel of dimension N with additive Gaussian noise and average power constraint. The encoder 

32 will use a bandwidth expansion mapping in order to encode a; to a point s{x) G respecting the constraint 

33 E [s{x)s{xY] < P. The decoder will recover an estimate x for the sent value attempting to minimize the mean 

34 square error (mse) E[{X ~ X)'^]. The source will be assumed to have pdf with support [0, 1) and s([0, 1)) will be 

35 referred as the signal locus. If s is a continuous mapping, then the signal locus is a curve in M^. 



X 




Z = {Zi, . . .,Z2n) 

z, = N{0,a) 
Fig. 1: Communication system 



36 This communication model is carefully analyzed in ||6l Ch. 4]. Given a radius p > 0, let Ep denote the event 

37 ||z|| < /9. The mse can be split into two terms such that: 

mse = E[{X - X f\Ep]P{Ep) + E[{X - l)2|£;^]P(i?p. (1) 

38 For high values of the SNR P/cr^, P{Ep) w and the mse is well approximated by E[{X - XY\Ep] (the 

39 goodness of the approximation depends on the SNR and p). In this case, called the low noise regime, it is possible 

40 to show ||6l that: 

E[{X - X f] « / f{x) \\s{x)\\-^ dx E^^^[{X - X)% (2) 
Jo 



41 where f{x) is the pdf of the source. Most of the analysis in Section VI will be done under the low noise regime. 



42 If the source is uniform, then an optimal parametrization for s{x) yields 

2 

i?low[(^-^)'] = ^, (3) 

43 where L is the curve length. 

44 A geometric picture of the possible effects of the noise on the estimated value is given by Figure 1. Under the 

45 low noise regime, all the calculations can be essentially done by approximating s{x) for its tangent line and the 

46 system will be described by Equation (|2]l. However, this is not true for large errors. Decoding to the wrong fold of 
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47 the curve yields a poor performance of the system. This suggests that our objective is to choose long curves with 

48 a large distance between its folds. 

49 In our constructions, the signal locus will be a collection of curves in M.^ rather than a single curve (this can 

50 be viewed as a "piecewise continuous" curve). Our analysis is mainly focused on uniform sources, although the 

51 scheme can be easily adapted to any source with limited support and even to Gaussian sources, by means of a 

52 companding function ||6l. 



Fig. 2: Small errors and large errors. 



III. Background 



54 A. Flat Tori 

55 The unit sphere 5^^^^ C can be foliated with flat tori fl], jSl as follows. For each unit vector c 

56 (ci,C2,..,cjv) e ]R^,Q > 0, and u = (ui, uj, ■ • • , wjv) G K^, let $c : ^ be defined by 



$c (m) = ci COS — , sin — , . . . , Civ cos , sin . (4) 

\ \ Ci Ci J \ Cat cjv / / 

57 This A^-periodic function is a local isometry on its image, the torus Tc, a flat A^-dimensional surface contained 

58 in the unit sphere S'^^^^. Tc = $c(^^) is also the image of the hyperbox: 

Vc ■■= {u e M^; < < 27rcJ, l<i<N. (5) 

59 Note also that each vector of S*^^"^ belongs to one, and only one, of these flat tori if we also consider the 

60 degenerated cases where some Ci may vanish. Thus we say that the family of flat tori Tc and their degenerations, 

61 with c = (ci, C2, .., Cat), ||c|| = 1, > 0, defined above is a foliation on the unit sphere of S"^^"^ C M^^. 

62 It can be shown (Proposition 1 in 18.]) that the minimum distance between two points, one in each flat torus 

63 and Tc, is 

1/2 



d{Tc,Tt) = 11c - 6|1 = ij2{c, - b,fj . (6) 
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The distance between two points on the same torus Tc given by 

64 is bounded in terms of Hit — by the following proposition. 

65 Proposition 1. l\8^ Let c = (ci, C2, .., cat), ||c|| — 1, and let u,v £ Vc- Let A = ||u — v\\ and 6 

66 Then 

sin sin — 

-A < — ^A < ,5 < ^A < A 

67 where = min Ci. 

68 

69 B. Curves 

70 As pointed out earlier in f9l, some important properties of a curve from a communication point of view are 

71 its stretch and small-ball radius. Given a curve s : [a, h] i-> M^, the Voronoi region V{x) of a point s{x) is 

72 the set of all points in which are closer to s{x) than to any other point of the curve. If H{x) denotes the 

73 hyperplane orthogonal to the curve at s{x), then the maximal small-ball radius of s is the largest r > such that 

74 B.r{s{x)) n H{x) C V{x) for all x E [a, b], where Br{s{x)) is the Euclidean ball of radius r centered at s{x). This 

75 means that the "tube" of radius r placed along the curve does not intersect itself The stretch S{x) is the function 

76 |ls(a;)|| where s{x) is the derivative of s{x). If ||s|| does not depend on x (as it is the case of the curves considered 

77 here) we will refer to the stretch as simply S. The length of a curve is given by ^^S{x)dx. In this paper we will 

78 be interested in curves with large length and small-ball radius. 

79 IV. The Torus Layer Scheme 

80 The design of those curves in this section is essentially divided in two parts. First, we choose a collection of M tori 

81 on the surface of the Euclidean sphere S"^^'^ at least 25 apart. The approach for doing this is via discrete spherical 

82 codes in M^. Second, we show a systematic way of constructing curves on each layer, via projection-lattices in 

83 M^^^. Finally, we give a description of the whole signal locus and summarize the encoding process. 

84 For this torus layer scheme, the unit interval will be partitioned into M intervals of different length, and each 

85 of them mapped into a curve on one of the layers. It is worth noticing that, for the special case M ~ 1, if we 

86 choose the torus associated to the vector c = ■ ■ ■ A) to encode the information, then the scheme proposed 

87 here is exactly the one analysed in ||9j. As we shown next, for M > 1, the curves presented here outperform the 

88 ones presented in f9\ (see also fT^) in terms of length and small-ball radius. 

89 A. Torus Layers 

90 Given a target small ball radius 6, the first step of our approach is to define a collection T = T2, • • • , Tm} 

91 of flat tori on S'^^^^ such that the minimum distance ^ between any two of them is greater than 26. This step is 



^ WMu) - Mv)\\. 

(7) 
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92 equivalent to design a TV-dimensional spherical code with minimum distance 25 and such that its M codewords 

93 have non-negative coordinates. 

94 Each point c e SC+ defines a hyperbox Vc ^ and hence a flat torus Tc in the unit sphere S"^^^^. By (|6| we 

95 can assert that the distance between two of those tori is at least 25. 

96 There are several ways of constructing spherical codes that can be employed here, e.g. ISj, H or even on layers 

97 of flat tori as introduced in JS). 

98 We present next an example of a family of discrete spherical codes in that will be useful in Section 



VI 



99 Example 1. Let 

100 be a non-negative N -dimensional unit vector. 

The set SC^^t) ~ {(T(c(t)) : a G 5* at} of all permutations of c{t) defines a N -dimensional spherical code with 
non-negative coordinates, minimum distance equals to 

d{t) = 



101 and cardinality M = N\. This set can be used in the first step of the construction described in this section. 

102 Note that 

limc(t) = ^(l,...,l):=e, (9) 



103 and, for each t > 0, the codewords of SCc(t) o'"^ equidistant from e, which is associated to the maximum volume 

104 torus Tg. 

105 Moreover, for any fixed dimension N and do CzR sufficiently small, there exists t > such that d{t) = do. This 

106 follows by observing that the equation 

do = 



has a positive root for 



0<do< 



y/{N -1)N{2N -1)' 

108 Although the construction above is not the best one in terms of number of points for fixed d, it has the advantage 

109 of being highly symmetric and having a closed formula for the spherical code minimum distance. 

110 B. Curves on Each Torus 

111 Let Tc be a torus represented by the vector c G SC+ as defined in the previous section. On the surface of Tc we 

112 will consider curves of the form: 

St.(x) - $c(a;2™), (10) 
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113 where C = diag(ci, . . . , cjv), u = uC = (ciUi, . . . , cjvujv), is given by ^ and x e [0, 1]. 

114 Provided that u £ Z^,gcd(Mi) = 1, those curves are closed (a (wi, . . . ,MAr)-type knot), and due to periodicity 

115 and local isometry properties of $c their lengths are 27r They are also the image through $c of the intersection 

116 between the set of lines l-nW , W — {ux + n . n = nC and n E Z^}, and the box Vc- For c = e, these are exactly 

117 the curves analyzed in ||9l- 

118 Let rc{u) be the minimum distance between two different lines in W. We have: 

rc{u) := min min \\ux — {ux + n)\\ 

= min min||ua; — n| 

(11) 

n^ku,k£Z 

= min \\Pi,±{n)\\ , 

119 where P^±{n) denotes the orthogonal projection of n onto the hyperplane which is given by the standard 

120 projection formula 



P,.(n)=n^/.--j. (12) 

121 If we consider ~ c-{L ® . . . cjvZ, the rectangular lattice generated by matrix C, then Vciu) is the length 

122 of shortest non-zero vector of the projectiorj^of Ac onto . Due to (|7]l, the small-ball radius 8u,c of St^ can be 

123 bounded in terms of Tc(u) as follows: 

2c,sin('^)<^„,<2sin('^), (13) 

124 where = minQ and ci > 0. Thus, for small values of S^^e, we have Su,e ~ 7rrc(w). Our goal is to choose u in 

125 order to maximize rdu). In addition, we also want to reach a contrary objective, which is the one of maximizing 

126 the arc length l^.e = 27r ||it|| of 8t^. 

127 The proportion (density) of the volume of the "tube" of radius TTrdu) inside Vc is precisely the packing density 

128 of the lattice P^±{Ac), the projection of Ac onto u-^, which is given by |5|: 




A(P„. (Ac)) = V^._i ( ^];_[^nJ^ j < A^.„i (14) 

129 where Ajy-i is the density of the best lattice in dimension {N — 1) and Vn^i is the volume of the {N — 1)- 

130 dimensional unit sphere. For the case when all entries of c are equal, Ac is equivalent to and it was shown in 

131 Q that we can make the above bound as tight as we want. We will show in Section V that this is also true for an 

'in general, the projection of a lattice A onto a subspace H is not a lattice unless certain special conditions are met, e.g., when H-^ is 
spanned by primitive vectors of A f5l. This will be always the case in this paper, since H = ii^ for a primitive vector u. 
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132 arbitrary c i.e., that projections of the rectangular lattice can also yield dense lattice packings and therefore we 

133 can construct curves on the flat torus with the parameters arbitrary close to this bound. 

134 Example 2. Let N — 2. Consider the local isometry 

^ I \ ( Ui . Ui U2 . U2\ 

3>c(w) = ci cos — , ci sm — , C2 cos — , C2 sm — (15) 

V Ci Ci C2 C2J 

135 on the flat torus Tc and the line segment given by v{x) = xv ^ x(27rMiCi, 27ru2C2), ,Ui,M2 G Z = x2ttu and 

136 < a; < 1. The curve St^{x) will be the composition ^{v{x)) and we have 

rdv) \\v\\ = 27rciC2 ^ rdv) = —^==== (16) 



2^2 



137 This curve in will turn around ui-times the circle obtained by its projection on the flrst two coordinates, 

138 whereas turning around U2-times the circle of radius C2 given by its last two coordinates (a type (ui,M2) knot on 

139 the flat torus Tc). In this case, it is easy to calculate rc{u) since the density of the projection lattice is equal to 

140 one. In Figure^it is illustrated the curve (^1,1/2) — (4,5), with ci > C2. 

141 C. Encoding 



142 



Let T = {Ti,...,Tm} be a collection of M tori as designed in Section IV-A For each one of these tori, 

43 let St^{x) — ^{x2TTUk), {k = 1,2,..., A/) be the curve on Tk, determined by the vector ilk ( flOl l and consider 

44 L = X^jli ^i' where Lk is the length of st^- 

Now split the unit interval [0, 1] into M pieces according to the length of each curve: 



and consider the bijective mapping 



[0,1) = /iU/2,...U/m, where 

for A: = 1 M. 



fk : h ^ [0, 1) 

— hjl — ■ 

145 Then the full encoding map s can be defined by 

s(a;) := STfc(/fe(a;)), if a; e 4. (17) 

146 The stretch of s will be constant and equal its total length L and the small-ball radius of s is the minimum 

147 small-ball radius 5 of the curves st^., provided that the distance between any pair of torus in T is at least 25. 

148 To encode a value x within [0, 1] we apply the map ( fTT] !. The signal locus will be a set of M closed curves, 

149 each one lying on a torus layer Tk and defined by a vector Uk. This whole process is illustrated in Figure [3] 

150 If the source is uniformly distributed over [0, 1], the encoding process presented above is a proper one, since all 

151 subintervals will be equally stretched (see SectionjVl]). For other applications, however, it could be worth considering 

152 another partition. 





Fig. 3: Encoding Process 



153 D. Decoding 

154 Given a received vector y e M?^ , the maximum likelihood decoding is finding x such that: 



i = arg min llj/ — s(x)||. 

a;e[0,l] 

155 Since exactly solving this problem is computationally expensive we focus on a suboptimal decoder 



56 For 7^ 7i = \/y2i-i + J/li' '^^ write 



157 where. 



vi y2\ ^ ( y2N-i J/2JV 

= 71 — 7iV , 

71 7i / V 7JV 7iv 

^1 . ^1 \ / On . On 

y = 71 cos — , sm — , . . . , 7JV cos — , sm — 

7l 7l / \ JN JN 



arccos ( 1 < i < N. 



It 

158 The process of finding the closest layer involves a iV-dimensional spherical decoding of 7 = (71,72, ■ ' ' ,7iv) 

159 which has complexity 0{MN). 

Let Ci = (cii, Ci2, • • ■ , CiAf) be the closest point in €E SC+ to 7 and 

- f f di . Oi\ f On . (^N 

J/i = ca cos — , sm — , . . . , c^at cos — , sm — 
V V 71 71/ V 7^ 7JV 

be the projection of y in the torus Tc-, i.e., 

\\y-y^\ < \\y-v\\,yv&T,^. 
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160 From now on, we proceed the process by using a slight modification of the torus decoding algorithm |j9l applied to 

161 the A^-dimensional hyperbox Va- The complexity of this algorithm is given by 0{N ||iti||j^), where Ui is the vector 

162 that determines the curve sxi- Hence, if M = 0(maxi ||u,;||]^), the overall complexity of the process described in 

163 this section will be 0(A^max,j the same as for the torus decoding. 



164 V. A Scaled Lifting Construction 

165 A. The construction 

166 The Lifting Construction was proposed in |j7l as a solution to the problem of finding dense lattices which are 

167 equivalent to orthogonal projections of along integer vectors ("fat-strut" problem). In this section we adapt 

168 that strategy in order to construct projections of the lattice Ac which approximate any {N — 1) dimensional lattice 

169 (hence the densest one) with the objective of finding curves in approaching the bound ([14]). We adopt here the 

170 lattice terms as in ||2l. For our purposes, the proximity measure for lattices will be the distance between their Gram 

171 matrices, as in Q. This notion measures how close a lattice is to another one up to congruence transformations 

172 (rotations or reflections). 

We consider the dual of a lattice A € M.^ 

A* = {x e span(A) : {x,y) e Z Vy e A}, 

173 where span(A) is the subspace spanned by a basis of A. Now let Ac — c\L © . . . ® catZ be the rectangular lattice 

174 generated by the diagonal matrix C (and with Gram matrix CC* = C^). By scaling Ac, we may assume that ci = 1. 

175 With this condition, if P„i (Ac)* denotes the dual of the projection of Ac through a vector u = (1, U2C2 . . . ,unCn) 

176 (ui e Z), then a generator matrix for P^j_(Ac)* is given by 



M = 



f-U2 1/C2 ... \ 

-«3 I/C3 ... 



(18) 



\-Un ... 1/C„.y 

177 what can be derived as a consequence of Prop. 1.3.4 f?]. In what follows we derive a general construction of 

178 projections such that P^± is arbitrarily close to a target lattice A e M^^^. 

179 Theorem 1. Let Ac — Z ® C2Z. . . catZ, q G M. Let A be a target lattice in M.^^^ and consider a lower 
ISO triangular generator matrix L* — (l*^) for A*. If A*^, w G N the sequence of lattices generated by the matrices 



T* 



C2 
C2 



Cn-1 



\ 







(19) 



then: 
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182 (i) i*„ = P^± (A)* for some ueR^ and 

183 (ii) (1/u>2)L;l;* ^ i*L** as w^oo. 

184 /n other words, for large w, A will be approximated (in the sense of the Gram matrices) by projections of Ac- 

185 Proof Through applying elementary (integer) operations on L*^^ we can put it on form ([T8| for some integers 

186 U2,.--,Un depending on w, hence LJ^ is a generator matrix for P^±{A)*, proving the first statement. For the 

187 second statement, we clearly have (l/w)L^ — > [L* 0] as w — > oo, where is the (rt — 1) x n all-zero column 

188 vector Therefore, {l/w^)LlL*J [L* 0] [L* 0]* = L*L*\ ■ 

189 Since the density of a lattice is a continuous function of the entries of its Gram matrix, it follows that the sequence 

190 of curves produced by the vectors u approaches the bound ( [T4| l as — > oo. 

191 Example 3. Consider the hexagonal lattice [2], which is the best packing in two dimensions and is equivalent to 

192 its dual. One of its generator matrix is 




193 We apply the construction above and reduce, through elementary operations, the matrix L*^ in order to put it on 

194 form ( |18| l. After re-scaling the rectangular lattice Ac, we find the sequence of vectors 

Uw = {ci,-2wc2,{2wlwV3c2/ci\ ~ 10)03). 

195 The projections of c\L © C2Z © C3Z onto will be, up to equivalence, arbitrarily close to A2 when w 00. 



196 B. Comparisons: Curves in 

197 Here we compare our approach of constructing curves on torus layers with the curves constructed in Q and Q 

198 in terms of length for given small-ball radii. Given (5 > 0, we first consider a set of flat tori associated to a spherical 



199 code in M.^, with minimum distance greater than 26, as described in Section IV-A Through the first inequality of 

200 ( [T3| l, for each torus, we can find rc in order to guarantee that the curves on the flat tori will have small-ball radius 

201 at least 5 (this is also done in the case of the curves on the torus c = e used in |9|). We then look for the larger 

202 element of the sequence of vectors that produces a projection with minimum distance at least rc- 

203 In Figure |4] the first curve from the bottom to the top represents the exponential sequence ID . The second one 

204 is obtained by directly applying the Lifting Construction [Jj to the hexagonal lattice and the last one displays the 

205 total length associated to our scheme. 



206 

207 
208 



VI. Mean Square Error Analysis 

The first step of the encoding process described in Section IV-C| involves a choice of partition of the interval 
[0, 1]. If the source to be transmitted is uniformly distributed over [0, 1], it is intuitive that the best choice in terms 
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209 of mse is the one such that the encoding map will have constant stretch. This fact is formalized in the next theorem 

210 and holds for uniform sources. 

Theorem 2. Let /i, . . . , Im be a partition of the interval [0, 1). Let {si, . . . , Sm} be a collection of non-overlapping 
curves on S"^^^^ of length Li, . . . , Lm and with domain Ii, . . . , Im respectively. Suppose the encoding rule 

s{x) = Sj{x) if X <E Ij 

211 is used. Then, under the low noise regime, the minimum e* of Ei^^[(X — X)"^] over all partitions {/i, . . . , /„j} is 



212 achieved by choosing Ij such that 



= (20) 



213 and satisfies e* — jL^. 



Proof: First note that 



M 



M 

and if the noise is small, Eig^[{X - X)'^\X e Ij] = a'^\Ij\'^ / L]. Therefore 
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M I r 13 / I 7- I \ 3 



(a) 



j=i J j=i ^ ' ' ^ 



(T 



j=i ^ / 

215 where the inequahty (a) is due to the strict convexity of the function g{x) = x'^ for a; > 0. Equality in (a) only 

216 holds when \Ij\/Lj — \Ik\/ Lk^k,]. This condition, together with X]j=i ^ 1' implies Equation (|20|. 

217 ■ 

218 It was proved in |9| that for fixed N and cr^, the mse decays with order 0{1/P^) provided that the parameter 

219 a [9 Sec. VI. A] is chosen suitably (i.e., the mse decays non-linearly with the power). In what follows we prove 

220 that the lengths of the curves constructed from the torus layer scheme are about 0{N\) greater than the scheme 

221 proposed in |9l, even if we use the improvements provided by the Lifting Construction Q. This will yield an 

222 asymptotic behavior comparable to the 0{1/P^) in |9|, but with a better performance when N increases. 

223 Let N be fixed and let c{t) be given by Equation (|8]l. Recall that e = l/\/iV(l, . . . , 1). For p > sufficiently 

224 small, there is t such that the tori in the set SCc(t) have minimum distance 2p. Moreover, there is li e ^c(t) 

225 arbitrarily close to the bound ([14]) i.e., such that the curve in each torus has length: 



Lc(t) = MM = ^-]^n W " ^0 

i—1 \ 1 / 



226 The same vector u can be used for any c £ SCc{t) just by interchanging its coordinates, yielding the same 

227 parameters for all tori. 

Since we have iV! tori in the set SCc{t), the total length Ltl produced is given by: 

228 For the Lifting Construction and the same small-ball radius, we have 

Llc = L, = 2n\\u\\ = - s, 

229 As p ^ 0, we can make both ei and 62 vanish, and then 

230 Now, by making t sufficiently small (but keeping the small-ball radius greater than p), we have 



^ A^! as t ^ 0. 



231 In other words. 
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232 where the approximation is good when we allow the small-ball radius to be sufficiently small (this is the case when 

233 the SNR is high). 

234 VII. Simulations 

235 We now present some simulation results in order to illustrate the performance of our proposed scheme. In all 

236 simulations we consider a source uniformly distributed over the interval [—1/2, 1/2]. 

237 As a benchmark we chose the spherical code scheme (V&C) presented in ||9l, since this scheme outperforms 

238 the one in [10] for uniform sources, as shown in their simulations. We also plotted the curves corresponding to a 

239 simple linear modulation. 

240 The simulation results corresponds to 50000 samples with 1 to 6 expansion rate {N = 3) and the performance 

241 is measured in terms of 1/mse as a function of the SNR = P/a^. 

242 In Figure |5] it is displayed the torus layer scheme performance for i\f = 6 tori designed as in Example 1 with 

243 t — 0.6, and projection vector u = (1,2, 198). The spherical code scheme V&C was simulated for a — 18. The 

244 parameters were chosen such that both schemes reach the asymptote for the same SNR (i.e., have approximately 

245 the same small-ball radius). As expected, our scheme outperforms (V&C) in the low noise regime (about 13dB 

246 higher in this example). 



120 




% 25 30 35 40 45 

SNR[dB] 

Fig. 5: Both schemes reach their asymptote for the same SNR. The torus layer scheme outperforms (V&C) in the 
low noise regime. 
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^% 22 24 26 28 30 32 34 36 38 40 
SNR [dB] per channel 



Fig. 6: Both schemes have about the same performance in low noise regime. Torus layers reaches the asymptote 
for a smaller SNR. 

247 Figure |6] shows the simulations results for M = 6 tori, also designed as in Example 1 with t = 1.75, and 

248 projection vector u = (1,4, 34). The spherical code scheme V&C was simulated also for a = 18. These parameters 

249 were chosen such that both schemes achieve about the same performance in low noise regime (same asymptote). 

250 In this case, our scheme reaches the asymptote for a smaller SNR (5dB in this example). In both graphics it is 

251 possible to see that the linear modulation is a good choice only for very small SNR. 

252 We finish this section with a remark on the implementation. Since the curves constructed from our method 

253 are closed, if the sent value is one extreme value of the interval Ij = [Lj^i/ L, Lj / L] it may happen that the 

254 decoder wrongly interpret it as the other extreme. In other words, if Lj^i/L + e is sent, then with a non-negligible 

255 probability, the decoder will output Lj/L. To prevent these errors, the encoder needs to "shrink" the interval Ij 
25$ by a factor a E (0, 1), such that the curve is "opened", separating the encoded extremes values apart. However, 

257 this process deteriorates the low-noise performance, so that a needs to be calibrated. In our examples, it sufficed 

258 to take any a between 0.7 and 0.8. This issue was briefly commented in |[9l- 

259 VIII. Conclusion 

260 The problem of transmitting a continuous alphabet source over an AWGN channel was considered through an 

261 approach based on curves designed in layers of flat tori on the surface of a (2A^)-dimensional Euclidean sphere. 

262 This approach explores connections with constructions of spherical codes and is related to the problem of finding 

263 dense projections of the lattice ciZ ® . . . catZ. 

264 This work is a generalization of both the scheme proposed in f9\ and the Lifting Construction in fT^. As a 

265 consequence, our scheme compares favorably to previous works in terms of the tradeoff between total length and 

266 small-ball radius, which is a proper figure of merit for this communication system. 

267 In spite of the improvements in terms of length versus small-ball radius, the constructiveness, homogeneity and 

268 overall complexity of the decoding algorithm are features preserved from |9|. Simulations also show a good scaling 
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269 of the mse with the channel SNR. Further questions include applications to sources that do not have limited support, 

270 or when it is not possible to employ a companding function (e.g., Gaussian mixtures), and generalization to other 

271 bandwidth expansion ratios then 1 : 2^^^ (the general case is to consider K to N expansion mappings, when 

272 KK < N). 
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